Quantum opening of the Coulomb gap in two dimensions 
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For a constant density of spinless fermions in strongly disordered two dimensional clusters, the 
energy level spacing between the ground state and the first excitation is studied for increasing 
system sizes. The average indicates a smooth opening of the gap when the Coulomb energy to 
Fermi energy ratio r s increases from to 3, while the distribution exhibits a sharp Poisson-Wigner- 
like transition at r s « 1. The results are related to the transition from Mott to Efros-Shklovskii 
hopping conductivity recently observed at a similar ratio r s . 
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In disordered insulators, a crossover M in the tem- 
perature dependence of the resistivity p(T) is induced 
by Coulomb interactions from the Mott variable range 
hopping law (p(T) — pu expfTo/T) 1 / 3 in dimension 
d = 2) to the Efros-Shklovskii behavior (p(T) = 
Pes expfTEs/T 1 ) 1 ^ 2 )- The long range nature of the in- 
teractions leads to a dip in the single particle density 
of states, and the assumption that single electron hop- 
ping dominates the transport leads to this change in the 
resistivity. However, a single electron hop may reorga- 
nize the location of the other particles, inducing complex 
many particle excitations. This makes the Coulomb gap 
problem difficult and gives us the motivation to study 
the first quantum excitation above the ground state. For 
d = 2, the strength of Coulomb interactions is very often 
given in units of the Fermi energy by the dimensionless 
ratio r s . For strong disorder, the first excitation energy 
is expected to become larger when r s increases, and to 
decay as 1/L (instead of l/L 2 for free electrons) when 
the system size L increases. But one cannot estimate 
the threshold rf where the Coulomb gap opens without 
taking into account all complex of many body quantum 
processes. Considering spinless fermions in 2d strongly 
disordered clusters, we confirm from numerical calcula- 
tions that the gap opens at a value r^ w 1.2, and we point 
out that this is indeed for a similar value (r s w 1.7) that 
a change in the hopping conductivity from Mott hopping 
to Coulomb gap behavior has been reported for an 
electron gas created at a GaAs/AlGaAs heterostructure. 
For a statistical ensemble of clusters, we have calculated 
the many body states at the mean field level given by 
the Hartree-Fock approximation and we have added the 
effects of the residual interaction. Keeping constant the 
carrier density n e , and increasing the size L, the average 
gap between the ground state and the first excitation be- 
haves as 1/L a , with a decreasing from 2 to 1 when r s 
increases from to 3. Another remarkable effect of the 
interaction is to yield a sharp transition for the gap distri- 
bution: it tends to Poisson or to Wigner-like distributions 
for small or large r s respectively at the thermodynamic 



limit. A critical threshold r^ ss 1.2 is characterized by 
a scale invariant gap distribution, reminiscent of the one 
particle problem j3| at a mobility edge. However, it is 
only the distribution of the first spacing which exhibits 
such a transition, the distributions of the next spacings 
remain Poissonian and are essentially unchanged when 
r s varies. Eventually, we discuss the implications for the 
hopping conductivity and we confirm that the transition 
for the gap takes place at a smaller r s than rf « 4 — 5 
where a change in the topology of the persistent currents 
carried by the ground state has been observed . 

We consider a disordered square lattice with M = L? 
sites occupied by N spinless fermions. The Hamiltonian 
reads 
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where cj (c$) creates (destroys) an electron in the site 
i, the hopping term t between nearest neighbours char- 
acterizes the kinetic energy, m the site potentials taken 
at random inside the interval [— W/2, +W/2], m = c\a 
is the occupation number at site i and U measures the 
strength of the Coulomb repulsion. The boundary con- 
ditions are periodic and is the inter-particle distance 
for a 2d torus. If a* B = h 2 e/(m,*e 2 ), to*, e, a and 
n s = N/(aL 2 ) denote respectively the effective Bohr ra- 
dius, the effective mass, the dielectric constant, the lat- 
tice spacing and the carrier density, the factor r s is given 
by: 
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since in our units Ti /(2m*a 2 ) 
n e = N/L 2 . 
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U and 



In this study, a large disorder to hopping ratio W/t — 
15 is imposed for having Anderson localization and Pois- 
sonian spectral statistics for the one particle levels at 
r s = when L > 8. We study N = 4,9 and 16 particles 
inside clusters of size L = 8, 12 and 16 respectively. This 
corresponds to a constant low carrier density n e = 1/16. 
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A numerical study via exact diagonalization techniques 
for sparse matrices is possible only for small systems j|, 
and does not allow us to vary L for a constant density. 
We are obliged to look for an approximate solution of 
the problem, using the Hartree-Fock (HF) orbitals, and 
to control the validity of the approximations. One starts 
from the HF Hamiltonian where the two-body part is 
reduced to an effective single particle Hamiltonian ||-||] 



(3) 



where (...) stands for the expectation value with respect 
to the HF ground state, which has to be determined self- 
consistently. For large values of the interaction and large 
system sizes the single-particle problem (||) is still non- 
trivial, since the self-consistent iteration can be trapped 
in metastable states. This limits our study to small r s 
and forbids us to study by this method charge crystal- 
lization discussed in Q at a larger rY ~ 12 . 

The mean field HF results can be improved using a 
method PJlO[| known as the configuration interaction 
method (CIM) in quantum chemistry jy}. Once a com- 
plete orthonormal basis of HF orbitals has been calcu- 
lated {Hhf^o) = ea|"0a) with a — 1,2, . . . , L 2 ), it is 
possible to build up a Slater determinants' basis for the 
many-body problem which can be truncated to the Nh 
first Slater determinants, ordered by increasing energies. 
The two-body Hamiltonian can be written as 



Hnt = \ Yl QZp d l4> dsd 7> 



(4) 



look if the results vary when Nh increases. In the worst 
case considered (L — 16, r s — 2.8) the accuracy in the 
first four spacings can be estimated of the order 5% when 
N H = 2x 10 3 . 



Therefore the CIM method allows to study low en- 
ergy level statistics for r s < 3. However, its accuracy is 
not sufficient to determine a small change of the ground 
state energy when the boundary conditions are twisted 
(i.e. the persistent currents). 
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FIG. 1. CIM result (full circle) compared to HF approx- 
imation (empty circle) for N = 9 and L — 12. Mean gap 
< Aq > and gap distribution at r s = 4.5 (insert). 



with 
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s a/3 



U 



(5) 



and S a = V'aOOcjjO). One gets the residual interac- 
tion subtracting Eq. [| from Eq. [|. This keeps the two- 
body nature of the Coulomb interaction, and if N 2 
it is still possible to take advantage of the sparsity of the 
matrix and to diagonalize it via the Lanczos algorithm. 

We have first compared HF and CIM results. Labelling 
the levels by increasing energy and studying an ensem- 
ble of 10 4 samples, we have studied the first spacing 
Ao = Ei — Eq . The role of the residual interaction can 
be seen in Fig. [I]. When r s > 1, the residual interaction 
reduces the mean gap, and slightly changes the distribu- 
tion. The CIM results agree with the results given from 
exact diagonalization with an accuracy of the order 2% 
when one takes into account the — 10 3 first Slater 
determinants when r s = 5 and L — 8. This means that 
a basis spanning only 0.2% of the total Hilbert space is 
sufficient for studying the first excitations. For larger L, 
exact diagonalization is no longer possible, but one can 
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FIG. 2. Size dependence of the average gap (first 
spacing < Ao >cx L - "' 7 " 3 '). From bottom to top: 
r 3 = 0, 0.6, 0.9, 1.1, 1.4, 1.7, 2.3, 2.8. Insert: a(r s ) (circle, char- 
acterizing < Ao > and square, average over i = 1 — 3, char- 
acterizing < Aj > ). 
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We have calculated the first energy levels for different 
sizes L. The first average spacing < Ao > calculated 
for an ensemble of 10 4 samples is given in Fig. ||. It 
exhibits a power law decay as L increases, with an expo- 
nent a given in the insert. One finds for the first spacing 
that a linearly decreases from d = 2 to 1 when r s in- 
creases from to 3. This proves a gradual opening of 
the mean Coulomb gap. The next mean spacings depend 
more weakly on r s , as shown in Fig. |^. 
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FIG. 3. Gap distribution P(s) for r s — (empty circle) 
and r s = 2.8 (full circle) when L — 16, compared to Pp(s) 
and Pw(s). Insert: size invariant P(s) at ~ 1.2 ; L — 8 
(circle), 12 (square) and 16 (diamond). 

For r s = 0, the distribution of the first spacing s = 
Ao/ < A > becomes more and more close to the Pois- 
son distribution Pp(s) = exp(— s) when L increases, as 
it should be for an Anderson insulator. For a larger r s , 
the distribution seems to become close to the Wigner 
surmise P\v{s) — (tts/2) exp(— 7rs 2 /4) characteristic of 
level repulsion in random matrix theory, as shown for 
r s = 2.8 and L = 16 for instance. To study how this P(s) 
goes from Poisson to a Wigner-like distribution when r s 
increases, we have calculated a parameter r/ which de- 
creases from 1 to when P(s) goes from Poisson to 
Wigner: 



var(P(s)) — var(JWs)) 
var(Pp(s)) — var(Pv^(s)) ' 



(6) 



where var(P(s)) denotes the variance of P(s), 
var(Pp(s)) = 1 and vax(P w (s)) = 0.273. In Fig||, 
one can see that three curves r](r s ) characterizing the 
first spacing for L = 8, 12, 16 intersect at a critical value 
w 1.2. For r s < the distribution tends to Poisson 
in the thermodynamic limit, while for r s > it tends to 
a Wigner-like behavior [Q. At the threshold rf, there 
is a size-independent intermediate distribution shown in 



the insert of Fig. ||, exhibiting level repulsion at small s 
followed by a exp(— as) decay at large s with a « 1.52. 
This Poisson- Wigner transition characterizes only the 
first spacing, the distributions of the next spacings be- 
ing quite different. The insert of Fig. || docs not show 
an intersection for the parameter r\ calculated with the 
second spacing. 




FIG. 4. Parameter r/(r s ) corresponding to the first spacing 
Ao at L — 8 (circle), 12 (square) and 16 (diamond). Insert: 
rj(r s ) for the second spacing Ai. 

The second excitation is less localized than the first one 
when r s = 0, since the one particle localization length 
weakly increases with energy. This is only for L = 16 
that the distribution of the second spacing becomes close 
to Poisson without interaction, and a weak level repul- 
sion occurs as r s increases. The observed transition, and 
the difference between the first spacing and the follow- 
ing ones is mainly an effect of the HF mean field. For 
the first spacing, the curves rj calculated with the HF 
data are qualitatively the same. At the mean field level 
the first excitation is a particle-hole excitation starting 
from the ground state and requires an energy of the or- 
der U I L, with fluctuations around this mean value. The 
second excitation is again a particle-hole excitation start- 
ing from the ground state. The energy spacing between 
the first and the second excited state is given by the dif- 
ference of two uncorrelated particle-hole excitations and 
a Poissonian distribution follows naturally. For r s > 7 
the Gaussian-like HF distributions for Ao become more 
Wigner-like when the residual interaction is included. We 
point out that in metallic quantum dots also, the first ex- 
citation is statistically different from the others, as shown 
by numerical studies [p^[ within the HF approximation. 

The existence of a critical r s value for the opening of 
the Coulomb gap can be understood similarly to p3|. 
The single particle density of states around the Fermi 
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energy E F is given by p(E) « |2? — E F \/U 2 and the 
gap size A g = \E g — £7f| can be estimated from the 
condition p(E g ) w p, with p w 1/W mean density of 
states for W 3> t, obtaining A g » U 2 /W. According 
to Fermi golden rule, the inverse lifetime of a Slater de- 
terminant built from electrons localized at given sites is 
r t » t 2 (l/W0(A7£ 2 ), with N/{WL 2 ) density of states 
directly coupled by the hopping term of the Hamilto- 
nian (Q). Therefore at zero temperature quantum fluc- 
tuactions melt the Coulomb gap for T t ~ A s , giving 
r s w m 1. We conclude that a crossover from Efros- 
Shklovskii to Mott hopping conductivity is expected not 
only increasing temperature but also increasing carrier 
density, as observed in 

To measure possible derealization effects, we have cal- 
culated the number of occupied sites per particle £ s = 
N/^2iP 2 where pi — (^ol^il^o) is the charge density 
of the ground state at the site i. Around r s w 1.2 
and after ensemble average, the maximum increase of 
£ s compared to r s — is negligibly small (2%). These 
are mainly the distribution and the average value of 
the first excitation energy which exhibit noticeable ef- 
fects. This matters for the hopping conductivity. The 
usual argument is to consider the length L(T) where 
exp[— (2L/£(r s ) + Ao(r s )/fcT)] is maximum with the lo- 
calization length £(r s ) w £(Q). If one takes for Ao(r s ) its 
average value w (A + Br s )/L a{ - r ^ (see Fig. |TJ> , one ob- 
tains for the hopping resistivity a smooth and continuous 
crossover from Mott to Efros-Shklovskii hopping, given 
by: 



where 



p(T) oc exp 



T(r s 



T(r s ) 
T 
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Br., 



(7) 



(8) 



This prediction neglects the sharp transition in the dis- 
tribution of A at , which could be better included by 
considering a more typical value for Ag(r s ) than its av- 
erage, for instance obtained from the value sj, for which 
J Q b p(s)ds = b, with 6=1/2 for instance. This will in- 
troduce a sharp discontinuity at in T(r s ). 

In summary, we have analyzed the Coulomb gap statis- 
tics for spinless fermions in a strongly disordered squared 
lattice when r s < 3. On one hand, we have found a 
sharp interaction- induced transition at rf « 1.2, char- 
acterized by a scale invariant distribution. Around the 
critical point, the gap distribution tends to Poisson or to 
Wigner-like distributions respectively at the thermody- 
namic limit. This effect is present at the HF mean field 



level, the residual interaction weakly shifts rf and im- 
proves the Wigner-like character of one of the limits. On 
the other hand, the exponent a characterizing the aver- 
age gap smoothly decays from 2 to 1 in this range of r s 
values. The average gap is substantially reduced from its 
HF value by the residual interaction. We associate this 
transition to a crossover in the hopping resistivity inside 
an insulating phase. 
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